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ABSTRACT 


Data from space from infrared imaging systems and SIRS are 
already providing information of great value in determining sea sur- 
face temperature and the atmospheric structure in cloudless areas 
over the oceans. Although some interpretations differ, it is clear that 
simultaneous measurements of radar sea return and passive micro- 
wave temperature will provide estimates of the wind speed, and per- 
haps wind direction, over the oceans, especially in cloudless areas, 
for a wide range of wind speeds. This report is concerned with the 
problem of integrating the data that would be obtained by a spacecraft, 
especially one with a combination radar-radiometer, into global analy- 
sis procedures for meteorological, wave and oceanographic pre- 
dictions. 

A computer based procedure that defines the mean vector wind 
as a function of elevation from the sea surface to the gradient wind 
level using an adequate theory of the structure of the planetary boundary 
layer has been developed. This procedure employs conventional ship 
reports, both those transmitted synoptically and those collected clima- 
tologically. Examples of the different fields that are produced are 
given. 

A global wave forecasting procedure has been developed. Im- 
proved meteorological predictions obtained from the system to be de- 
scribed below would provide high accuracy wave forecasts in the two to 
four day time frame. Examples of the wave specifications produced 
for the North Pacific are shown. Data from space can in principle 
verify the wave heights given by this model as of a given time of ob- 
servation. 
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So as to study areas of poor conventional data coverage, the 
wind analyses have been used to generate data on wind speeds such as 
would be obtained from a spacecraft. A small fraction of the avail- 
able ship reports for the North Pacific Ocean have been used to simu- 
late the data density for the South Pacific Ocean. A combined know- 
ledge of low center positions from cloud patterns, coastal and island 
surface pressures, a few scattered ship reports, and the simulated 
spacecraft data permit the recovery of the vector wind field over the 
oceans and the surface atmospheric pressure patterns. The surface 
pressure pattern is in turn an important parameter in the interpre- 
tation of SIRS data. A combined analysis of the surface pressure field 
and the SIRS data would yield information for the entire atmospheric 
structure in the initial value problem for a global weather forecast. 

The conditions in the interior of the oceans, that is, currents, 
temperatures, and salinities, are controlled by what happens at the 
coasts and at the air-sea interface. Fluxes of radiation, heat, mo- 
mentum and water across the air-sea interface and river runoff and 
ground water discharge at the coasts can be used as inputs to an ocean 
circulation model. Presently available data from space and the data 
that could be obtained by future spacecraft should make it possible both 
to specify and predict conditions inside the oceans in much the same 
way as the weather is predicted at present. Efforts to develop such 
systems are described. 


PRESENTLY AVAILABLE DATA FROM SPACE 


Quantities of importance to both meteorology and oceanography 
are presently being sensed from space. Examples are SIRS data, 
infrared sea surface temperature data, and clouds and cloud patterns. 
These measurements provide atmospheric soundings, sea surface 
temperature fields, in which temperature fronts in the oceanic struc- 
ture can be detected, and information on the location of low centers 
and fronts. 

McClain (1970) has described how sea surface temperature 
fields can be obtained for four to five day average conditions over 
both hemispheres, using presently available infrared imagery. The 
motion of the clouds reveals most of the ocean surface during this 
period of time, so that realistic fields can be constructed. The sea 
surface temperature does not vary rapidly so that four or five day 
averages are permissible for some applications. 

The ATS and NOAA imagery of clouds, such as shown in the ATS 
photograph in Fig. 1, show how fronts and low centers can be located. 
The surface low pressure centers are undoubtedly located a little to 
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the east and equatorward of the tips of the spiraling cloud patterns 
shown in this image. Low centers over the ocean are often missed 
completely or located poorly on the basis of conventional data, and 
such images already clearly locate lows as accurately or more ac- 
curately, depending on their location, than ship reports. The long 
arching bands locate fronts. Cold fronts are often carried forward 
by extrapolation for many days over the ocean on the basis of very 
poor data in conventional analyses. 

SIRS data as described by Wark and Hilleary (1969) can be used 
to define the constant pressure surfaces for the multi-layer numerical 
weather prediction models currently in use, as for example described 
by Shuman and. Hovermale (1968). It is important to note that a better 
definition of the pressure at the surface of the ocean would make the 
SIRS data even more valuable because it would reduce the error in the 
specification of the heights of constant pressure surfaces for numerical 
prediction procedures. 


RADAR AND MICROWAVE RESULTS 


Results on radar sea return for five different frequencies and on 
passive microwave emission over a wide range of frequencies lend 
further support to the concept of using a combination, active radar 
passive microwave instrument to learn about the winds and waves over 
the oceans. Ross et al. (1970) have shown that passive microwave 
emission is sensitive to the presence of whitecaps and that the micro- 
wave temperature is related to the foam on the sea surface for high 
winds and waves. Hollinger (1970) in this third annual review, shows 
that passive microwave emission at low winds is sensitive both to the 
structure of the sea surface for high wave numbers and to whitecaps, 
with the whitecaps being of lesser importance than the sea surface 
structure for low winds. 

The report by Guinard (1970), also in this annual review, indi- 
cates that radar sea return is weakly dependent on winds (and waves) 
for high winds, and that 13.3 gigahertz seems to be a good frequency 
to study radar sea return. The results of Guinard are not in agree- 
ment with the results of Moore and Pierson (1970) for 13.3 gigahertz, 
where a stronger dependence of radar sea return on wind at higher 
winds was found, based on a revised analysis of the available radar 
sea return data from several NASA missions. It is important to note 
that recent work by Valenzuela et al. (1970), using some of the same 
radar sea return data reported by Guinard (1970), can be interpreted 
to indicate a wind speed dependence even for higher winds. Pierson 
and Moore (1970) have discussed these matters In considerable detail. 
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The report by Krishen (1970), shows that the problems of cali- 
brating and correcting the 13.3 gigahertz scatterometer of NASA 
MSC have been solved and that the Mission 119 data and the data 
from previous missions will form a consistent set of o°(0) curves 
for a wide variety of wind speeds that wTll be good to ± ~ db. Some 
data for Mission 119 has already been prepared and released. A 
preliminary analysis of these data shows that the value of (r°(0) 
agree with previously obtained data, as described by Moore and 
Pierson (1969), for both the upwind -downwind condition and for 
the crosswind condition. 

One important point about radar sea return is that it is direction 
dependent with all workers reporting a substantial difference between 
upwind -downwind and crosswind measurements. In contrast, the 
passive microwave data do not appear to be sensitive to direction, 
either the direction of the waves on the sea surface or the direction 
of the wind that generated the waves. For cloudless conditions, it 
may be possible to combine sea surface temperature measurements 
with the infrared, passive microwave emission measurements and 
radar sea return measurements to obtain both information on the 
wind speed and partial information on the wind direction. 

Radar sea return varies from a maximum in either the upwind 
or downwind direction to a minimum at the crosswind direction. If 
the wind speed could be determined from the passive microwave 
measurements, then if the radar is looking in the upwind or down- 
wind direction, the radar sea return will be a maximum and the wind 
direction will be known as either the direction \ or the direction 
X + 180°. Similarly if the radar sea return is a minimum for the 
given velocity, then the wind direction is again known with an indeter- 
minacy of 180°. If the radar sea return value lies between the maxi- 
mum for upwind -downwind and the minimum for crosswind, this im- 
plies some angle, say <p , between the upwind -downwind direction of 
the wind and the radar beam. In turn, this implies that there are 
four possible wind directions that would be given by x ± </> and 
X ± (f> + 180°. The overall pattern of winds around low centers and a 
few scattered ship reports should then provide a high quality analy- 
sis of the vector wind and the pressures in the planetary boundary 
layer. 


For other conditions such as viewing through stratus clouds 
with a high liquid water content or rain, the passive microwave sig- 
nal will be needed to calibrate the radar signal, but it should still be 
possible to obtain data on the wind speed over the sea surface under 
these conditions. 
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COMPUTER BASED PROCEDURES FOR THE USE OF 
DATA FROM SPACECRAFT 


Once the winds and the pressures in the planetary boundary 
layer have been determined, they provide the information needed for 
the planetary boundary layer in numerical weather prediction models. 

If the pressure field can be determined accurately at the sea surface, 
it provides a basic input to the SIRS data that then makes it possible 
to describe the elevations of the constant pressure surfaces through- 
out the entire troposphere as described by Smith and Fritz (1969). 

Thus a combination of SIRS data and data from an active radar, pas- 
sive microwave instrument could almost by itself, along with infrared 
sea surface temperature data, define the initial value problem for a 
numerical weather prediction for the entire southern hemisphere, given 
a few scattered ship reports and cloud imagery. For the northern 
hemisphere, the spacecraft data will improve the specification of the 
initial value problem compared to presently available data. This im- 
plies improved forecasts of the weather and improved wave forecasts. 

It is, however, necessary to develop the computer procedures 
for utilizing the data that would become available from a combination 
passive microwave active radar instrument on a spacecraft in a polar 
orbit. This has been a major effort of our group during the past few 
years and has resulted in new techniques for analyzing the wind and 
pressure fields in the planetary boundary layer by means of conventional 
data. This work has also resulted in the development of a highly supe- 
rior numerical wave specification and wave forecasting procedure and 
in the development of a system that simulates the data one might get 
from a radar radiometer and shows how it can be incorporated into a 
conventional analysis so as to provide an analysis of the planetary 
boundary layer in areas of low conventional data coverage. Finally, 
all of these developments point toward the possibility of computing what 
occurs and will occur in the interior of the oceans, just as is now done 
for the atmosphere so as to obtain true oceanographic predictions. 


ANALYSIS OF THE PLANETARY BOUNDARY LAYER 


If the meteorological and oceanographic data are not obtained in 
real time, the data density is substantially increased because many 
ships send in their logs of the weather a month or so late to the National 
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Weather Records Center. This data base is thus quite complete for 
most of the Northern Hemisphere, with the amount of ships available, 
being approximately twice the number that are available on the day- 
to-day real time teletype services. With this augmented data base, 
Cardone (1969) has developed a computer based procedure for defin- 
ing the winds over the surface of the oceans. By means of this analy- 
sis, the winds as a function of height from the surface to the gradient 
wind level can be defined over the Northern Hemisphere oceans. The 
analysis produces the vector wind at twenty meters, the Monin-Obukov 
mixing length, the vector friction velocity at the sea surface, which 
in turn can easily be transformed to the vector wind stress, and the 
wind at the gradient wind level since the thermal wind is known. Part 
of the intermediary data needed consists of the air-sea temperature 
difference field. Examples of the various fields produced by this com- 
puter-based procedure are shown in Figs. 2 through 8. Fields such 
as these are available for the Northern Hemisphere every six hours 
for more than one year. The printouts are distorted and they should 
be mentally stretched so that the decimal points would lie on the sides 
of squares. The title on each figure describes it. The grid is four 
times the density of the conventional NWP analysis. The analysis has 
not been smoothed to the same extent as is usual for hemispherical 
analyses. 

The air-sea temperature field as shown in Fig. 6 is an im- 
portant field for describing the variation of the wind with height in 
the first meters above the sea surface. This field was obtained from 
conventional sources. Infrared imagery can supplement the present 
data sources for determining this field, and SIRS data could give a 
few data points along the subsatellite track for the air temperature 
near the sea surface. Above line 105, which is toward the equator, 
there is not enough information available to specify this field but 
the values are generally small and have been assumed to be zero. 

The field for the inverse Monin-Obukov stability length shown 
in Fig. 7 provides the required parameter for the computation of 
the wind profile over the water. If it is zero, the wind profile is 
logarithmic. Positive values mean a stable atmospheric stratification 
and a weaker wind stress on the surface for the same wind at ane- 
mometer height along with a different variation of the wind with height. 
Conversely negative values mean an unstable stratification and a 
greater stress on the sea surface for the same wind at anemometer 
height. 

Fig. 8 shows the magnitude of the friction velocity on the sea 
surface. The directionof the friction velocity is the same as in Fig. 

5 for the wind at 19.5 meters. This stress field, which is directly 
obtained from this field, is an important field for numerical models 
of the wind driven ocean circulation. 
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One of the major purposes of the work by Cardone was to produce 
wind fields that would yield more accurate wave specifications; that is, 
a computation given the observed winds for all past times of what the 
waves are at a particular time, and more accurate wave forecasts, 
that is, a computation of what the waves will be, given the winds and 
waves up to a certain time, and forecasts of the winds for a day or 
so into the future. The growth of the wave spectrum in a given fre- 
quency range depends on the wind speed, the atmospheric stability, the 
fetch, the duration of the wind, and on spectral components that have 
propagated into an area from other areas. Each of these factors varies 
from point to point over the ocean. All of these factors are properly 
included in the wave specification and wave forecasting procedures to 
be discussed next in this paper in terms of the theories of Phillips 
(1957), Miles (1957), and Phillips (1966). 

There is, however, a part of the wave spectrum that is satur- 
ated over a given frequency band. Further input from the wind in 
adding energy to the wave spectrum for these bands, only causes the 
production of whitecaps and breakers, so that this energy is event- 
ually dissipated in turbulence. Cardone (1969) developed a way to 
compute the amount of energy dissipated in this way and called this 
quantity, the whitecap production index. This whitecap production in- 
dex was then compared with the measurements of Monahan (1969) and 
with the measurements in the North Sea described by Ross et al. 

(1970). The theory for the computation of the whitecap production in- 
dex is shown in Fig. 9. The computation of this index is a straight- 
forward extension of present results in wave forecasting and specifi- 
cation. There should eventually be a strong interrelationship between 
the whitecap production index and the work of Austin (1970) as de- 
scribed at this Third Annual Earth Resources Review. 

Given the specifications and forecasting procedures to be de- 
scribed in the next section of this paper, and given the sea surface 
temperature fields as described by McClain, it should not be difficult 
to remove the effects of fetch, duration, atmospheric stability and 
sea surface temperature, so as to determine the wind speed simul- 
taneously with radar measurements. This would provide one of the 
quantities discussed above for obtaining information on the vector wind. 


WAVES 


Numerical wave specification and numerical wave prediction 
techniques have been developed. The report by Pierson, Tick and 
Baer (1966) discusses some of the theory that is used. Bunting (1970) 
has shown how well the predictions can be expected to verify. 
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Fig. 10 shows the land-sea table for the North Pacific and Fig. 11 
shows the land-sea table for the North Atlantic Ocean that will be 
used as a basis of printing out such quantities as wave height and 
wave direction for the hindcasts (and eventually forecasts) that will 
be produced with this computer program. 

The first successful run of this improved and highly sophisti- 
cated wave specification and wave forecasting system for the North 
Pacific was achieved during the week of December 7th, and the week- 
end of December 12 and 13, 1970. The results, as given in Fig. 12, 
show what the waves over the North Pacific would have been like if 
the sea surface had been flat calm and then the winds had started to 
act on the sea surface and had blown for six hours. After seven more 
days, or so, the waves over the North Pacific will be quite realistic 
since the observed winds are used to compute what the waves were 
like. 


The waves are not very high, and the values shown illustrate 
the problem of start-up time in numerical wave specification models. 
It is impossible to measure the two dimensional wave spectrum over 
the while North Pacific so as to obtain the equivalent of the con- 
ditions for a more classical initial value problem. The wave spectra 
specified by these procedures, after a week or so, given good wind 
fields, yield significant heights and frequency spectra that verify well 
against actually measured waves. 

It is quite possible that the data from the radar altimeter part 
of SI 93 on Skylab A will provide information on wave height that can 
be compared with the heights predicted by this model. The results 
of Yaplee et al. (1970) are quite encouraging in this connection as 
they lend strong support to the assumptions made by Pierson and 
Mehr (1970) in deriving the theoretical form for an average altimeter 
pulse being scattered back from the sea surface. 

These wave forecasting and wave specification procedures were 
developed for deep water. Waves near coasts are modified by the 
effects of refraction over irregular submarine topography. It is im- 
portant to note here the work of Chao (1970) in which the theory of 
wave refraction at a caustic was solved and of Chao and Pierson (1970) 
in which the theory was verified in a wave tank. It should now be 
possible to predict wave conditions at all important coasts as well as 
in deep water. 
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INTEGRATION OF THE DATA FROM A RADAR RADIOMETER 
INTO A PLANETARY BOUNDARY LAYER ANALYSIS 


It is also necessary to develop ways that the data from a radar 
radiometer can be integrated into conventional data so as to produce 
useful analyses of the wind and pressure fields in the planetary boundary 
layer. To study this problem, the data given by the work of Cardone 
has been assumed to provide an accurate description of the winds over 
the North Pacific ocean for a particular synoptic map time. Simulated 
passes of a polar orbiting spacecraft are made over the North Pacific, 
and it is assumed that the magnitude of the vector wind at 19.5 meters 
can be determined from the combined radar radiometer data on a 
grid of points corresponding to repeated cross track scans. This 
provides a dense field of wind speeds. It is further assumed that the 
location of low centers can be determined from space photographs of 
clouds as described above, that conventional coastal data and island 
data will be available and that a few scattered ship reports will also 
be available, along with a field of the air-sea temperature difference 
from, say, SIRS and infrared sea-surface temperature fields. The 
assumptions simulate what will be available from a combination of 
spacecraft mounted sensors and conventional data for the Southern 
Hemisphere, for example, over the South Pacific ocean. Only about 
twenty or thirty ships from the as many as three hundred to five 
hundred that actually report for the North Pacific are selected at 
random from the total list of ships to simulate the Southern Hemisphere 
conditions. 

The results of a computer based procedure developed at New York 
University by Mr. Leonard Druyan that uses these isolated ships, 
simulated radar radiometer data reduced to wind magnitudes alone as 
each grid point, island and coastal data, and cloud data are shown in 
Fig. 13. The recovery of the pressure field for the limited number of 
ships used is quite good. Since this figure was prepared added modi- 
fications of the computer procedure have led to even better results. 

It is planned, once the computer based procedures have been 
further improved, to compute the pressure fields for a year of data, 
once per day, so as to produce statistics on how well the system can 
be expected to operate. The data produced by Cardone, as described 
above, will also be used to obtain information on the relative frequency 
of different wind speeds over the Northern Hemisphere oceans, for a 
further evaluation of the potential accuracies of the proposed radar 
radiometer system. 

Also, the concept that some information on the wind direction 
can be inferred from the space data has yet to be evaluated in terms 
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of the use of limited ship data and spacecraft data in the Southern 
Hemisphere. Limitation of the possible wind directions to either two 
or four directions at a given point from combined radar radiometer 
data, can provide enough information when combined with infor- 
mation on the global circulation to define the pressure pattern even 
more accurately. The simulation described above that uses only wind 
speed may actually be an overrestriction on the potential system. 

The root mean square error of the pressure field in Fig. 13 is al- 
ready low enough to make the pressure field obtained by this method 
useful in applying SIRS data to the determination of the heights of 
constant pressure surfaces over the ocean. 


PREDICTION OF THE OCEAN CIRCULATION 


The time may soon be upon us when it will be necessary to pre- 
dict conditions inside the ocean, just as weather conditions are now 
predicted for the atmosphere. The quantities to be predicted at a 
minimum would be the current speed and direction and the tempera- 
ture and salinity at a number of levels inside the ocean. The con- 
ditions in the interior of the ocean are however, controlled by effects 
at the ocean atmosphere interface and by runoff at the coasts. Models 
of the oceans, such as those of Bryan (1969), show clearly that if the 
fluxes of heat, water, radiation and momentum across the air-sea 
boundary could be defined over the oceans, and if river runoff could 
be computed, the motions and conditions in the interior of the ocean 
would respond realistically to these forcing functions. It would not 
be necessary to measure the currents, temperatures and salinities 
in the interior of the ocean in order to predict what goes on in the 
interior of the ocean, if the conditions at the boundaries are correctly 
specified throughout past time. 

Work toward this goal is also underway. Ways to use space 
data to compute the solar radiation absorbed at the sea surface are 
under development and are being correlated with XERB-1 data. The 
XERB-1 buoy is an experimental buoy located off our east coast that 
is a prototype of the b uoys that may be deployed by the National Data 
Buoy System of NOAA. There may be many of these buoys deployed 
in the Northern Hemisphere. They could be the surface truth sites 
for an eventual combination spacecraft-buoy system. Ways to 
compute the infrared loss from the sea surface have been developed. 
The SIRS soundings would be used to define the atmospheric structure 
over the oceanic points for which the infrared loss would be com- 
puted. However, the data density may not be great enough since the 
method works only along the subsatellite track. The stress fields 
are known from the work of Cardone, as described above. Other 
work should make it possible to apply real time forcing boundary 
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conditions to the finite difference equations of motion for the oceans 
as given by Friedrich (1970) and Okubo (1970), that will cause a real- 
istic ocean circulation throughout the ocean interior. The motions 
and conditions in the ocean will then vary in response to the actual 
weather as it occurred over the ocean surface. 


CONCLUSIONS 


It seems, therefore, that data from space could play an import- 
ant role in producing weather forecasts for the Southern Hemisphere 
of at least the quality now available for the Northern Hemisphere, in 
producing improved weather forecasts for the Northern Hemisphere 
bu augmenting the data base for the initial value problem in the North- 
ern Hemisphere, in producing forecasts of waves on the global ocean, 
and in producing forecasts of what is occurring, and will occur in the 
ocean interior. Some of the necessary computer based procedures 
for accomplishing these goals have been developed and others are 
now being developed in our present program. 
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Figure 1. ATS Image showing fronts and low centers near spiral cloud 
bands over the North Atlantic and South Atlantic Oceans. 


Asterisks represent land. 

Date: 12/1/66 Time: 0000 
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Figure 



Wind direction in degree* for the portion of the expanded (125 x 125) NWP grid system covering the western NortH Pacific Ocean. 
Asterisks represent land. 
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Figure 


Asterisks represent land. 
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Figure 9. The computation of the whitecap production index. 
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Fig. 13. Top : Pressure field analysis based on all available ships. Bottom: Pressure 
field analysis based on 17 ships and simulated wind speed information. 


